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A simple analytical model for dark matter halo structure and 

adiabatic contraction 
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A simple analytical model for describing inner parts of dark matter halo is considered. It is assumed that 
dark matter density is power-law. The model deals with dark matter distribution function in phase space of 
adiabatic invariants (radial action and angular momentum). Two variants are considered for the angular part 
of the distribution function: narrow and broad distribution. The model allows to describe explicitly the process 
of adiabatic contraction of halo due to change of gravitational potential caused by condensation of baryonic 
matter in the centre. The modification of dark matter density in the centre is calculated, and is it shown that the 
standard algorithm of adiabatic contraction calculation overestimates the compressed halo density, especially in 
the case of strong radial anisotropy. 
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1 Introduction 

The structure and evolution of dark matter halos is widely 
discussed in the recent years (e.g. pQ) A large number of 
papers has been devoted to the investigation of spatial 
density profiles of dark matter halos, p(r). both numer- 
ical [21 El and analytical 0]. Recently several papers 
appeared that described general properties of solutions 
of Jeans equation [SJ|S]. The Jeans equation for a speri- 
cally symmetric self-gravitating system reads as follows: 



d . 2 , 2/3 . 2 . GM{r) 



0. 



(1) 



are the radial and tangential velocity 

(7 2 • 

— „ * a is the Binney's anisotropy pa- 



dr 

Here ay 2 and a 2 
dispersion, j3 = 1 

rameter (systems with < f3 < 1 are radially anisotropic, 
with /3 < - tangentially anisotropic, and f3 = is the 
isotropic case). 

An additional assumption is the following empirical 
property of dark matter halos: the generalized phase- 
space- like density of dark matter particles in a halo is a 
power-law in radius: 

p(r)/o- r e (r)ocr-«, (2) 

which is established both in N-body calculations and 
some analytical models 6 . 

For a closed set of equations one also needs to spec- 
ify f3. Here two different approaches may be adopted: in 
[H] isotropy was assumed (/? = 0), and in jS] the follow- 
ing relationship between (3 and the logarithmic slope of 
density profile 7 = d In p(r) / dr was taken: 



This kind of relation has been proposed in [H] based on 
empirical analysis of different numerically treated situa- 
tions. 

These investigations demonstrated that under these 
assumptions there exist only few possibilities for a self- 
consistent halo model. The most realistic is the set 
of models with density profiles having inner and outer 
power-law asymptotes (so-called Zhao a(3j models 9 ). 

However, this approach, being very fruitful, does not 
describe explicitly the distribution function. It operates 
only with its second moment with respect to velocity, i.e., 
the velocity dispersion. The most complete description 
of halo structure can be made in terms of phase-space 
distribution function f(r,v). One attempt to deal with 
full velocity distribution function (VDF) was made in 
where an Eddington inversion formula was applied 
to obtain VDF for power-law density profiles. However, 
this study was limited to isotropic case. 

An approach based on certain physical assumptions 
regarding the form of the distribution function seems 
more promising. It avoids postulating empirical assump- 
tions like (j2H3J), which, nevertheless, still hold in many 
cases. It is convenient to express the distribution func- 
tion in terms of adiabatic invariants during slow evolu- 
tion of the gravitational potential. This allows to de- 
scribe easily the process of adiabatic contraction of dark 
matter halo caused by condensation of baryonic mat- 
ter in its centre, which occurs during the formation of a 
galaxy. This approach is adopted in the present paper. 

The paper is organized as follows. In the second sec- 
tion we describe two models of inner halo structure and 



13 = fa + 6(7 - 7o ) . 



(3) 
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explain physical arguments for the underlying assump- 
tions concerning the dark matter distribution function. 
In the third section we calculate the velocity anisotropy 
coefficient. Furthermore, in the fourth section we con- 
sider the halo response to the adiabatic compression and 
compare it with other studies. Finally, the conclusions 
are presented. 

2 The phase-space structure of model 

HALO 

Dark matter halos are extended complex objects which 
formed in different physical processes (growth of initial 
perturbations, collapse, hierarchical merging, and later - 
baryonic contraction) . We are intererested in the central 
area of a halo, which presumably is less influenced by 
merging and more - by condensation of baryons. Most 
of existing dark halo models assume power-law behaviour 
of density profile in the centre of a halo. 



p d (r) = Kr" )! , 1 < 7 < 2. 



(4) 



We shall adopt this assumption, having in mind that our 
analysis applies only to the inner area of a halo, interior 
to the radius where the logarithmic density slope starts 
to deviate from the constant value 7. The anisotropy 
coefficient is also taken to be constant in this area (0 < 
(3 < 1 , according to analytical calculations and numerical 
simulations ^3)- Under these assumptions, / depends 
on E also in a power-law form J2| ■ 

We write distribution function in terms of action- 
angle variables. The action variables are the angular 
momentum to, its projection onto z axis m z , and radial 
action defined as I r = § v r (r) dr. The distribution 
function is quasi-stationary and hence does not depend 
on angle variables, and due to spherical symmetry, nei- 
ther it depends on m z . The radial action in power-law 
potential is well approximated by the following expres- 



sion: 



I r = AE^T 



f 



m c (E) 



(5) 



This relation, according to ^1], emerges as a result of dy- 
namical mixing in the centre of a collapsing dark matter 
halo. Initial density peak, which gives birth to the grav- 
itationally bound object, is taken to be triaxial ellipsoid. 
Gravitational collapse and subsequent phase mixing, as 
shown in this paper, lead to formation of spherically sym- 
metric structure, in which particles have angular momen- 
tum directly linked to their initial distance from centre, 
and, hence, with their radial action. The quantity 1$ in 
expr. JfjJ) is related to orbital eccentricity of particles: 
la = corresponds to pure radial orbits, and — > 00 - 
to circular. 

The distribution function in model A is expressed as 
follows: 



f(I r ,m)=f E 1 / 2 6(m 2 -l 2 I r 2 ). 



(7) 



In the model B the distribution function is given by 
expression 



f(E, to) = f E» m- 2B , M = \ - \-l 



(1-B). (8) 



where m c (E) is the angular momentum of circular orbit 
with given energy |UI| . 

We consider two variants of inner halo structure, which 
differ in the dependence of distribution function on an- 
gular momentum. 

In the model A the angular momentum to of a par- 
ticle is proportional its radial action: 



— let Ir 



(6) 



This is the simplest generalization of isotropic model for 
the case of arbitrary radial anisotropy [TzIHB) . 

Therefore, the two models considered represent, in 
some sense, two limiting cases - narrow and broad dis- 
tribution by angular momentum for fixed energy. In 
what follows, we express particle energy through vari- 
ables I r and to according to QJ, and rewrite the distri- 
bution function in action-angle variables. 

Notice that we do not consider possible change of 
distribution function in the subsequent evolution, hier- 
archical clustering and merging, which is not significant 
for the inner parts of a halo, as shown in analytical cal- 
culations ^B] and numerical simulations |17| . 

Due to scale invariance, the expression is essen- 
tially satisfied, with a = (1 — Qj-)e. For 7 < 2 and e > 2 
the power-law index a is slightly greater than 2. This 
is similar to ESIM (Extended Secondary Infall models) 
described in [Hj , in which the density has also power-law 
form at small radii. 

3 Velocity anisotropy 

As long as the distribution function is given in action- 
angle variables, one can compute velocity distribution 
function at each radius, and, consequently, obtain radial 
and tangential velocity dispersion values ay 2 and a 2 ■ 
Obviously, they should have also power-law form a 2 oc 

r 2-7_ 
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In the model A the anisotropy coefficient (3 is related 
to the density slope 7 via the additional parameter 1$ 
(see Fig.l). This relation can be approximated as 

k 7- l-0.27Zo7- (9) 

Hence, for a given Iq the relation between 7 and (3 has 
linear form J3J), though the coefficients differ from [H]. 
The upper bound on /3 equals 7 — 1, in agreement with 

HE! 

In the model B the anisotropy coefficient (3 is identi- 
cal to B in expr. © ^2]. Notice that it is possible to 
have (3 > 7— 1, if we recall that power-law density profile 
(0J breaks down at large radii. In this case the velocity 
dispersion is determined by density profile cutoff scale, 
and depends on radius as a 2 oc r 7_2/3 . Upper bound on 
j3 is then 7/2, according to I n the model A it is 

not possible, since it follows from expr. (0 HJJ that the 
particle velocity at each radius lies in a limited range, 
proportional to ^/<5(r), and hence, for sufficiently small 
radii the cutoff scale does not determine the velocity dis- 
persion. 

4 Response to adiabatic contraction 

It is well-known that central parts of galaxies are domi- 
nated by baryonic matter. The response of dark matter 
halo to the baryonic infall usually is calculated using 
the adiabatic approximation. The common algorithm 
for computing compression is that of Blumenthal et al. 

in which dark matter particle's angular momentum 
is conserved. It is strictly applicable only to particles on 
circular orbits. In this method, the final radius rf of a 
dark matter particle is related to its initial radius r$ by 
the following expression: 

(Mf in , DM (rf) + M finMr (rf)) r f = M in . DM (ri) n , 

(10) 

where Mfi n- DM + ^fin,bar is the total mass of dark mat- 
ter and baryons after contraction, and Mi n ^oM is the 
initial dark matter mass inside a given radius. 

In general case of non-circular orbits, the radial ac- 
tion is also conserved, being an adiabatic invariant. How- 
ever, usually this does not help much, as the radial action 
is not known in general case. An improved method pro- 
posed in [20] involves the combination M{r) r, where r is 
the orbit-averaged radius of a particle. This yields less 
compressed halos, in better agreement with numerical 
simulations. Alternatively, an iterative method can be 
used for computing density profile and radial action in 
baryon-induced potential 21 . 




Fi gure 1 : Anisotropy (3 vs. density slope 7 relation in model A: 
dashed lines — initial anisotropy for different values of Iq ; 
solid lines - anisotropy after compression. Corresponding lines for 
the same lg cross at 7 = T = 1.8. Horizontal axis is the slope 
7 for original (uncompressed) halo and 7' for compressed, which 
depends on 7 via 1111 . Corresponding original values of 7 are 
shown in the bottom auxiliary axis. 
Thick line — upper limit j3 = 7 — 1; 
dotted line - Hansen & Moore relation 8 

In our simple models we are able to explicitly consider 
the adiabatic contraction, since we have our distribution 
function already written in terms of adiabatic invariants 
(provided that the characteristic timescale of the com- 
pression is much greater than dynamical timescale). All 
we need is to calculate radial action in new potential 
and recompute density profile. In many cases we may 
assume the baryonic density profile to be also power-law 
in radius: pb{r) = Kbt~ t , with index T < 2 (close to 
isothermal profile). This again applies only to central 
area of a galaxy, the bulge. We assume T > 7. 

It is quite obvious that the modified dark matter den- 
sity profile is also power-law: p' d {r) — K' r~ 7 . The new 
power-law index 7' is expressed as follows: 



The power-law index is essentially the same as in Blu- 
mcnthal's method, but the normalization coefficient is 
different. We introduce the quantity 77 = K' /K' std , the 
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Figure 2: r; — the ratio of our model's density normalization 
to that of the standard adiabatic contraction model, for different 
values of anisotropy coefficient /3. 

Dashed lines — model A, /3 = and 0.25 (from top to bottom); 
solid and dotted lines — models A and B correspondingly, for j3 = 0, 
0.25 and 0.5 (top to bottom), in model A the anisotropy coefficient 
/3' after contraction is used. 



ratio of our model's normalization to that of Blumen- 
thal's method (see Fig. 2). It appears that 77 < 1, in 
agreement with the results of [213] and [23? where it was 
found that halos with radial motion are compressed less 
than predicted by Blumenthal et al. algorithm. The 
quantity 77 decreases with decrease in 7 and increase in 
(3. Numerical simulations in [221 also confirm this re- 
sult; these authors explain it by the fact that for fixed 
[3 a greater value of 7 increases the dominance of circu- 
lar orbits over radial ones, and this makes Blumenthal's 
method more applicable. (Note that when j3 tends to 
7 — 1 in model A or to 7/2 in model B, the density nor- 
malization is more and more affected by density cutoff 
radius, so that values of 77 are underestimated). 

For quantitative calculations we take Milky Way bulge 
and dark halo as an example. The bulge density profile 
and normalization are taken from 123] and are the fol- 
lowing: p b = 0.6- 10 9 M Q /kpc 3 (r/lkpc)- r , T = 1.8. 
The dark matter halo is normalized to have density 0.3 
GeV/cm 2 at r=8.5 kpc |23J. (Note that it yields dif- 
ferent initial mass of dark matter inside 1 kpc, ranging 
from 4 • 10 8 M© for 7 = 1 to 3 • 10 9 M Q for 7 = 1.8. The 
baryonic mass within 1 kpc is 6 • 10 9 M©). 
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Figure 3: Ratio of compressed to initial dark matter mass within 
1 kpc: 

long dashed — compressed by standard (Blumenthal) algorithm; 
solid and dotted lines - models A and B correspondingly, for /3 = 
0,0.25 and 0.5 (from top to bottom); in model A the anisotropy 
coefficient /3' after contraction is used. 



We calculate the ratio of new to old dark matter mass 
within radius 1 kpc and plot it in Fig. 3. Again we see 
that with increasing (3 (more radial orbits), the actual 
compression ratio becomes smaller. The greater the dif- 
ference between 7 and T, the more is the increase in dark 
matter mass, i.e. the compression is more noticeable for 
halos with lower power-law index of initial density cusp. 

And finally, we calculate the velocity anisotropy pa- 
rameter f3' for modified halo (see Fig.l). In model A it 
appears to be generally greater than for corresponding 
uncompressed halo, meaning that velocity becomes even 
more radially biased. At the same time, the relation be- 
tween (3 1 and 7' becomes more flat (the same effect can 
be seen from Fig. 2 of jHj when comparing pure cosmo- 
logical simulations and dark matter with cooling baryons 
simulations). 

Having in mind that setting T = 7 = 7' effectively 
does not change profile, we approximate velocity aniso- 
tropy parameter as 



« /3(r) + 0.4(7' - r) . 



(12) 



In model B the anisotropy coefficient slightly decreases 
(at most by 0.04 for 7 = 1, /? = 0). We note that due to 
increase in /?' for model A it resembles more the model B 
in the reaction for adiabatic compression, if we plot the 
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compression coefficient as a function of j3', not (3 (solid 
and dashed lines in Figs. 2 and 3). 

5 Conclusion 

We have considered two simple analytical models for 
structure of central part of dark matter halo. In both 
models, the distribution function depends on energy in a 
power-law form, and the dependence on angular momen- 
tum is 5-function in model A and power-law in model B. 
The simplicity of the models allows us to calculate some 
interesting properties of the halo, namely: the velocity 
anisotropy and response to adiabatic contraction due to 
baryonic infall. These calculations help to get qualita- 
tive insight into halo properties and lead to the following 
conclusions: 

1. The more radially biased is the velocity, the less 
compressed is the halo. 

2. The shallower is initial density profile (less value 
of 7), the less is the degree of compression com- 
pared to that of BlumcnthaPs method. Standard 
model of adiabatic contraction systematically over- 
estimates the effect of contraction. 

3. The model predicts a moderate enhancement (2 to 
4 times) of dark matter mass in the bulge. The 
shallower is the initial profile, the greater the in- 
crease of dark matter mass. This result signifi- 
cantly increases the estimates of possible annihila- 
tion radiation flux from the center of the Galaxy 
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